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Abstract 

The reductive perturbation method has been employed to derive 
the Korteweg-de Vries (KdV) equation for smah but finite amphtude 
electron- acoustic waves. The Lagrangian of the time fractional KdV 
equation is used in similar form to the Lagrangian of the regular KdV 
equation. The variation of the functional of this Lagrangian leads to 
the Euler-Lagrange equation that leads to the time fractional KdV 
equation. The Riemann-Liouvulle definition of the fractional deriva- 
tive is used to describe the time fractional operator in the fractional 
KdV equation. The variational-iteration method given by He is used 
to solve the derived time fractional KdV equation. The calculations 
of the solution with initial condition sech(cx)^ are carried out. The 
result of the present investigation may be applicable to some plasma 
environments, such as the Earth's magnetotail region . 
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1 Introduction 



Because most classical processes observed in the physical world are non- 
conservative, it is important to be able to apply the power of variational 
methods to such cases. A method used a Lagrangian that leads to an Euler- 
Lagrange equation that is, in some sense, equivalent to the desired equation 
of motion. Hamilton's equations are derived from the Lagrangian and are 
equivalent to the Euler-Lagrange equation. If a Lagrangian is constructed 
using noninteger-order derivatives, then the resulting equation of motion can 
be nonconservative. It was shown that such fractional derivatives in the 
Lagrangian describe nonconservative forces [1, 2]. Further study of the frac- 
tional Euler-Lagrange can be found in the work of Agrawal [3, 4], Balcanu and 
coworkers [5, 6] and Tarasov and Zaslavsky [7, 8]. During the last decades. 
Fractional Calculus has been applied to almost every field of science, en- 
gineering and mathematics. Some of the areas where Fractional Calculus 
has been applied include viscoelasticity and rheology, electrical engineering, 
electrochemistry, biology, biophysics and bioengineering, signal and image 
processing, mechanics, mechatronics, physics, and control theory [9]. 

On the other hand, electron acoustic waves (EAWs) have been observed 
in the laboratory when the plasma consisted of two species of electrons with 
different temperatures, referred to as hot and cold electrons [10], or in an 
electron ion plasma with ions hotter than electrons [11]. Also its propagation 
plays an important role not only in laboratory but also in space plasma. For 
example. Bursts of broadband electrostatic noise (BEN) emissions have been 
observed in auroral and other regions of the magnetospherc, e.g. polar cusp, 
plasma sheet boundary layer (PSBL). see [12]. There are different methods 
to study nonlinear systems [13]. Washimi and Taniti [13] were the first to use 
reductive perturbation method to study the propagation of a slow modulation 
of a quasimonochromatic waves through plasma. And then the attention has 
been focused by many authors [14-15]. 

To the author's knowledge, the problem of time fractional KdV equation 
in coUisionless plasma has not been addressed in the literature before. So, our 
motive here is to study the effects of time fractional parameter on the electro- 
static structures for a system of unmagnetized coUisionless plasma consisting 
of a cold electron fluid and isothermal ions with two different temperatures 
obeying Boltzmann type distributions. We expect that the inclusion of time 
fractional parameter will change the properties as well as the regime of exis- 
tence of solitons. Several methods have been used to solve fractional differen- 
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tial equations such as: the Laplace transform method, the Fourier transform 
method, the iteration method and the operational method [16]. Recently, 
there are some papers deal with the existence and multiplicity of solution of 
nonlinear fractional differential equation by the use of techniques of nonlinear 
analysis [17-18]. In this paper, the resultant fractional KdV equation will be 
solved using a variational-iteration method (VIM) firstly used by He [19]. 

This paper is organized as follows: Section 2 is devoted to describe the for- 
mulation of the time-fractional KdV (FKdV) equation using the variational 
Euler-Lagrange method. In section 3, variational-iteration Method 9VIM) is 
discussed. The resultant time-FKdV equation is solved approximately using 
VIM. Section 5 contains the results of calculations and discussion of these 
results. 



2 Basic equations 

We consider a homogeneous system of unmagnetized coUisionless plasma con- 
sisted of a cold electron fluid and isothermal ions with two different temper- 
atures obeying Boltzmann type distributions. Such system is governed by the 
following normalized equations in one dimension [15]: 



d d 

—ne{x,t) + —[ne{x,t)ue{x,t)] = 0, (la) 
d d d 

[■Q^ + '^e{x,t)—]Ue{x,t)-—(t){x,t) = 0, (lb) 

x,t) -ne{x,t) +nii{x,t) +nih{x,t) = 0, (Ic) 



dx^ 

the two ions density nii{x, t) and nih{x, t) are given by: 

, , ,.— (f)(x, t), 

nii{x,t) = /xexp[ — J, (Id) 

nih{x,t) = i/exp[ — ^_^i^ll]. (le) 

In the earlier equations, ne{x, t) is the cold electron density normalized by 

equilibrium value Ueo , Ue{x, t) is the cold electron fluid velocity normalized by 
Ce// — (Teff/me)^, T^ff — ^Tu+^Ti ' temperature of low temperature 
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ion with initial normalized equilibrium density /i, is the temperature of 
high temperature ion with initial normalized equilibrium density i/, /3 = |f 
is the ions temperatures ratio, (j){x^ t) is the electric potential normalized by 
Tf,ff/e, is the mass of electron, e is the electron charge, x is the space 
co-ordinate normalized to the effective Debye length Xo^ff = (ijj^^f^)^ 
t is the time variable normalized by the inverse of the cold electron plasma 
frequency cu^^ , [cu^^ = ( "^^^"°° )^]- The neutrality condition reads fi + u — 1. 
Equations (la) and (lb) represent the inertia of cold electron and equation 
(Ic) is the Poisson's equation needs to make the self consistent. The two ion- 
densities are described by Boltzmann type distributions given by equations 
(Id) and (le). 



3 Nonlinear small-amplitude 

According to the general method of reductive perturbation theory [13], the 
slow stretched co-ordinates are introduced as: 

r^eh, ^ = e^x-Xt), (2) 

where e is a small dimensionless expansion parameter and A is the wave speed 
normalized by Ce//. All physical quantities appearing in (1) are expanded 
as power series in e about their equilibrium values as: 



ne(e,r) = 1 + em (e,r) + 6^712 (e,r) + e3n3(e,r) + ..., (3a) 
Me(e,r) = eMi(e,r) + e2u2(^,r) + eV3(e,r) + ..., (3b) 
0(e,r) = e0i(e,r) + e202(C,^) + e'V3(e,r) + -, (3c) 

with the boundary conditions that as |^| — )■ oo, = 1, = 0, = 0. 

Substituting (2) and (3) into (1) and equating the coefficients of like 
powers of e lead, from the lowest and second-order equations in e, to the 
following KdV equation for the first-order perturbed potential: 

^AMCrf-^ + B'-^^^O, (4a) 



where 
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In equation (4a), t) is a field variable, ,^ is a space coordinate in the 

propagation direction of the field and r G T(= [0, Tq]) is the time coordinate. 
The resultant KdV equation (4a) can be converted into time-fractional KdV 
equation as follows: 

Using a potential function V{^, r), where t) = V^{^, r) = $(^, r), 

gives the potential equation of the regular KdV equation (4a) in the form 

%(e, t)+A V^iC, t)v^^{C, t) + B %^5(e, r) = 0, (5) 

where the subscripts denote the partial differentiation of the function with 
respect to the parameter. The Lagrangian of this regular KdV equation (4a) 
can be defined using the semi-inverse method [20, 21] as follows: 

The functional of the potential equation (5) can be represented by 

J(y) = jd^j dT{y(e,r)[ciy^,(e,r)+C2AF^(e,r)^;^^(e,T)+C3Sy^^^^(e,r)]}, 

R T 

(6) 

where ci, C2 and C3 are constants to be determined. Integrating by parts and 
taking Vt-\r = V^Ir = V^\t = lead to 

J{V)^ j dij cir{V-(e,T)[-CiV-^(e,T)V-,(e,T)-^C2^V-/(e,T)+C3Sl-4(e,T)]}. 
R T 

(7) 

The unknown constants Q(i = 1, 2, 3) can be determined by taking the 
variation of the functional (7) to make it optimal. Taking the variation of this 
functional, integrating each term by parts and making the variation optimum 
give the following relation 

2ciy^.(e, r) + ?>C2AV^{i, r)%(e, r) + 2c,BV^^^^{i, r) = 0. (8) 
As this equation must be equal to (5), the unknown constants are given 

as 

ci = 1/2, C2 = 1/3 and C3 = 1/2. (9) 
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Therefore, the functional given by (7) gives the Lagrangian of the regular 
KdV equation as 



L(K, = -^^^(^,r)K(^,r) - lAVli^,r) + ^W/^(^,r). (10) 

Similar to this form, the Lagrangian of the time-fractional version of the 
KdV equation can be written in the form 



FioD^V, V^,V^^) = --[,D:V{C,T)]V^{^,r)--AVl{i,T) + -BVl^{C,T), 
< a<l, (11) 

where the fractional derivative is represented, using the left Riemann-Liouville 
fractional derivative definition as [16] 



k-1 < a<l,t e[a,b]. (12) 
The functional of the time-FKdV equation can be represented in the form 

J{V) - j dij drFioD^V, V^,V^^), (13) 

R T 

where the time-fractional Lagrangian F{qD^V, V^, V^^) is defined by (11). 

Following Agrawal's method [3, 4], the variation of functional (13) with 
respect to V{^, r) leads to 



SJ(V) = <lr{^S„D'^V + §.iV, + l^-^VJJ. (14) 

R T 

The formula for fractional integration by parts reads [3, 16] 



f dtf{t)aDtg{t) = f dtg{t)tDtf{t), f{t), g{t) e [a, h]. (15) 

J a J a 

where jD^, the right Riemann-Liouville fractional derivative, is defined by 
[16] 
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k-1 < a<l,t e[a,b]. (16) 
Integrating the right-hand side of (14) by parts using formula (15) leads 

to 



R 



dF d dF d'^ OF 



SJ{V) ^ jdil drWS^^) - -(— ) + —,{^)]SV, (17) 



where it is assumed that SV\t = SV\r = SV^\r = 0. 

Optimizing this variation of the functional J{V), i. e; SJ{V) = 0, gives 
the Euler-Lagrange equation for the time-FKdV equation in the form 

^„ , (9F , d ,dF^ d\dF ^ 

Substituting the Lagrangian of the time-FKdV equation (11) into this 
Euler-Lagrange formula (18) gives 



Substituting for the potential function, V^(^, t) = 0i(^, r) = $(^, r), 
gives the time-FKdV equation for the state function $(^,t) in the form 



lioD^^C, r) -r D^^ii, t)] + A $(e, r) $^ (C, r) + B ^^^^{i, r) = 0, (20) 

where the fractional derivatives o-D" and tD^q are, respectively the left and 
right Riemann-Liouville fractional derivatives and are defined by (12) and 
(16). 

The time-FKdV equation represented in (20) can be rewritten by the 
formula 

1 d 

- iD^m, t) + A $({, r)-$(C, r) + B r) = 0, (21) 
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where the fractional operator qD^ is called Riesz fractional derivative and 
can be represented by [4, 16] 



_ 1 1 d'' 
2r(A;-a)dt^7a 




k—a—l 



fir)], 



k-1 



< a <1, t e [a,b]. 



(22) 



The nonlinear fractional differential eqnations have been solved using dif- 
ferent techniques [16-20] . In this paper, a variational-iteration method (VIM) 
[21] has been used to solve the time-FKdV equation that is formulated using 
Euler- Lagrange variational technique. 

4 Variational-iteration method 

A general Lagrange multiplier method is constructed to solve non-linear prob- 
lems, which was first proposed to solve problems in quantum mechanics [21]. 
The VIM is a modification of this Lagrange multiplier method [22]. The basic 
features of the VIM are as follows. The solution of the hnear term of the prob- 
lem or the initial (boundary) condition of the nonlinear problem is used as 
initial approximation or trail function. A more highly precise approximation 
can be obtained using iteration correction functional. Variational-iteration 
method (VIM) [21] has been used successfully to solve different types of inte- 
ger nonhnear differential equations [22, 23]. Also, VIM is used to solve hnear 
and nonlinear fractional differential equations [24, 25]. This VIM has been 
used in this paper to solve the formulated time-FKdV equation. 

Considering a nonlinear partial differential equation consists of a linear 

part LU {x, t), nonlinear part NU {x, t) and a free term f{x, t) represented as 



where L is the linear operator and N is the nonlinear operator. According 
to the VIM. the (n+ l)th approximation solution of (23) can be given by the 
iteration correction functional as [24, 25] 



LU{x,t) + NU{x,t) = f{x,t), 



(23) 
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Un+i{x,t) ^Un{x,t)+ [ drX{r)[LUn{x,r)+NUn{x,r)-f{x,r)],n>0, 

Jo 

(24) 

where A(r) is a Lagrangian multiplier and Un{x,T) is considered 
stricted variation function, i. e; 6Un{x,T) = 0. Extreme the variation of 
the correction functional (24) leads to the Lagrangian multiplier A(r). The 
initial iteration can be used as the solution of the linear part of (23) or the 
initial value U{x,0). As n tends to infinity, the iteration leads to the exact 
solution of (23), i. e; 

U{x,t) = lim Un{x,t). (25) 

n— >oo 

For linear problems, the exact solution can be given using this method 
in only one step where its Lagrangian multiplier can be exactly identified. 

5 Time-fractional KdV equation solution 

The time-FKdV equation represented by (21) can be solved using the VIM 
by the iteration correction functional (24) as follows: 

Affecting from left by the fractional operator qD^"~^ on (21) leads to 



8 8^ 
- ^D^'--[A^^,r)-mr) + B ^$(e,r)], 

< a< 1, re [0,To], (26) 
where the following fractional derivative property is used [16] 



^D,-[^Dfm] = ^D^'m-t: ^D,^-=fm^a rli-a-'jy ^ < ^• 



(27) 

As a < 1, the Riesz fractional derivative ^-D^" ^ is considered as Riesz 
fractional integral that is defined by [16] 
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olrVit) = lioIrVit) + rhVit)] = ]^^j\r\t - f {t) , « > 0, 

(28) 

where Ql^°'f{t) and are the left and right Riemann-LiouvuUe frac- 

tional integrals, respectively [16]. 

The iterative correction functional of equation (26) is given as 

$„+l(e,T) = $„(^,r) + ^^dT'A(T'){^$n(e,T') 







la-2 



-0%' *n(e,r')k'=o^ 



T(q;-1) 

+ ^n(e, r')-^n{i. r') + £ ^$.(e, r')]}(29) 



where n > and the function $„(^,t) is considered as a restricted variation 

function, i. c; ^$„(,^,r) = 0. The extreme of the variation of (29) using the 
restricted variation function leads to 

5$„+i(e,T) = 5$„(e,T)+j^'dT'A(T')5^$n(e,T') 

This relation leads to the stationary conditions 1-|-A(r) = and ^A(r') = 
0, which leads to the Lagrangian multiplier as A(t') = —1. Therefore, the 
correction functional (29) is given by the form 



V{a-l) 

+ r')^$n(e, r') + B ^$.(e, r')]}(30) 
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where n > 0. 

In Physics, if r denotes the time-variable, the right Riemann-Liouville 
fractional derivative is interpreted as a future state of the process. For this 
reason, the right-derivative is usually neglected in applications, when the 
present state of the process does not depend on the results of the future 
development [3]. Therefore, the right-derivative is used equal to zero in the 
following calculations. 

The zero order correction of the solution can be taken as the initial value 
of the state variable, which is taken in this case as 

MC,r) = $(e,0) = Aosech\cC). (31) 
where ^40 = ^ and c—^ are constants. 

Substituting this zero order approximation into (30) and using the defi- 
nition of the fractional derivative (22) lead to the first order approximation 

as 



^i{^,t) = Aosech{c^y + 2Aoc smh{c^) sech{c^f 

*[4c'S + (^0^ - 12c'S) sech(cO'] ^, ^" (32) 

L[a + 1) 

Substituting this equation into (30), using the definition (22) and the 
Maple package lead to the second order approximation in the form 
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$2(C,r) = Aosech(c^)^ + 2Aoc sinh(cO sech(cO^ 
*[4c^B + {AoA - 12c^B) sech(cO^^ ^ 



r(a + l) 
+2Aoc^ sech(cO^ 

*[32c^52 + lQ>c^B{bAoA - 63c^S) sech(cO^ 
+2(3^^>1^ - nQ>Aoc^AB + leSOc^S^) sech(cO'' 

-7{AlA^ - A2Aoc^AB + 360c^£^) sech(cO" 



T(2q; + 1) 
+4AoC^sinh(cO sech(cO^ 
*[32c^B2 + 2Ac^B{AoA - Uc^B) sech(cO^ 
+A{AlA'^ - 32Aoc^AB + 240c^fi2) sech(c^)^ 
-5{AlA^ - 2AAqc^AB + lAAc^B^) sech(cO^] 

r(2a + 1) r^" 
*[r(a + l)]2r(3a + l)' 



(33) 



The higher order approximations can be calculated using the Maple or 
the Mathematica package to the appropriate order where the infinite approx- 
imation leads to the exact solution. 



6 Results and discussion 

Numerical studies have been made for a small amplitude electron-acoustic 
waves in an unmagnetized coUisionless plasma consisted of a cold electron 
fluid and isothermal ions with two different temperatures obeying Boltzmann 
type distributions. We have derived the Korteweg-dc Vries equation by using 
the reductive perturbation method [13]. The Riemann-LiouvuUe fractional 
derivative [16] is used to describe the time fractional operator in the FKdV 
equation. He's variational-iteration method [21] is used to solve the derived 
time-FKdV equation. 

However, since one of our motivations was to study effects of initial equi- 
librium density n of low temperature ion and time fractional order a on the 
existence of solitary waves. Our system can support two kinds of potential 
structure namely, compressive and rarefactive pulses. Depending on the sign 
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of the coefficient of the nonhnear term A, compressive sohton exists if A > 
while rarefactivc sohton exists if A < 0. 

In Fig (1), profiles of the bell-shaped rarefactive and compressive solitary 
pulses are obtained due to the change of the range of . Figure (2) shows that 
both the amplitude and the width of the compressive solitons increase with 
while both decrease for rarefactive solitons. Also, the time fractional order 
decreases the amplitude of the rarefactive and compressive solitons as shown 
in Fig (3). 

In summery, it has been found that amplitude and width of the electron- 
acoustic waves as well as parametric regime where the solitons can exist is 
sensitive to the initial equilibrium density of low temperature ion . Moreover, 
the time fractional order plays the role of higher order perturbation theory 
in changing the soliton amplitude. 

The application of our model might be particularly interesting in the 
new observations for the Earth's plasma sheet boundary layer region. We 
have stressed out that it is necessary to study the critical case for A — 0, the 
amplitude of the solitary pulse tends to infinity and the time-FKdV equation 
is not appropriate for describing the system. This is beyond the scope of the 
further work. 
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Figure Captions 

Fig (1): The electric potential $(,^,r) vs ^ and r for X — 1, v — 0.04, 
a = 0.5, (3 = 0.05, (a) = 0.2 and (b) = 0.3. 

Fig (2): The electric potential $(^,t) vs ^ and for A = 1, f = 0.04, 
a — 0.5, ^ = 0.05 and t — 5: (a) 3 dimensions and (b) 2 dimensions. 

Fig (3): The amplitude of the electric potential $(0,t) vs t and a 
for A = 1, v = 0.04, p = 0.05 and = 0.2: (a) 3 dimensions and (b) 2 
dimensions. 
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